Complete Two Loop Electroweak 
Contributions to the Muon Lifetime in the 

Standard Model 

M. Awramik 1 and M. Czakon 2 

1 Institute of Nuclear Physics, Radzikowskiego 152, PL-31342 Cracow, Poland 

2 Department of Field Theory and Particle Physics, Institute of Physics, 
University of Silesia, Uniwersytecka 4, PL-40007 Katowice, Poland 



Abstract 

An independent result for the two loop fermionic contributions to the muon life- 
time in the Standard Model is obtained. Deviations are found with respect to [1], 
which result in a shift of the W boson mass by approximately —1.3 MeV over the 
range of Higgs boson masses from 100 GeV to 1 TeV. Supplied with the bosonic 
contributions from [2,3,4], this shift, due to the complete electroweak contributions, 
varies from —2.4 MeV to —0.6 MeV. Additionally, a new test of the matching proce- 
dure defining the Fermi constant is presented, which uses fermion masses as infrared 
regulators. 



1 Introduction 

The muon lifetime is one of the key observables of today's particle physics. 
Not only is it measured very precisely, since the current experimental error is 
18 ppm [5], but can be described to competing accuracy within the Standard 
Model, giving rise to a strong correlation between the masses of the heavy 
gauge bosons. As a low energy process, the decay is expected to be governed by 
an effective interaction involving only the electron, muon and their respective 
neutrinos. The dynamics of the system should be corrected mostly by QED, 
whereas the electroweak interactions determine solely the size of the coupling 
constant. 

The history of the calculation of the electroweak corrections, in which we are 
interested here, is rather long. It started in the early eighties with the one 
loop contributions [6]. Subsequently, leading terms in the top quark 0(a 2 mf) 
[7] and Higgs boson 0(a 2 M]j) [8] masses were derived at the two loop level. 
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In the meantime, mixed electroweak and QCD corrections became available 
at order 0(aa s ) [9] and 0(aa 2 ) [10]. Recently, three loop leading terms in 
the top quark mass 0(a 3 m^) and 0(a 2 a s mf) have also been calculated [11]. 
As far as the pure two loop electroweak corrections are concerned, after it 
turned out that the subleading terms in the top quark mass expansion are 
comparable with the leading ones [12], complete fermionic and the Higgs boson 
mass dependence of the bosonic contributions have been evaluated [1]. The 
complete bosonic part has been done in [2,3,4]. It is the purpose of the present 
paper to present the result of a new independent calculation of the fermionic 
contributions and, after inclusion of the bosonic corrections, also of the full 
electroweak corrections. 

This work is organized as follows. In the next Section, we discuss the matching 
procedure and fermion masses as infrared regulators, which avoid ambiguities 
of the definition of gamma matrices in box diagrams in noninteger dimensions. 
Then, we present the results for the fermionic and full contributions and com- 
pare them with previous calculations by specifying the differences in the W 
boson mass prediction. Conclusions close the paper. 



2 Matching 



Due to a large number of very different mass scales, it is virtually impossi- 
ble to evaluate the muon decay lifetime directly within the Standard Model, 
without recourse to some approximation method. An elegant and systematic 
approximation is provided by the approach based on effective theories. The 
idea is to agree on some cutoff scale, below which all degrees of freedom are 
treated exactly, whereas the heavier fields are "integrated out" , which means 
that they generate effective interactions. It should not be surprising that one 
first discovers experimentally the effective theories, since the dependence on 
the heavier scales requires higher "resolution", i.e. higher energy. For precisely 
this reason, the effective theory governing muon decay, the Fermi Model, has 
been known much before the Standard Model. From this point of view, one 
should not consider that the Fermi Model is used in current calculations for 
historical reasons, but because it is the appropriate effective theory at this 
energy scale. 

The approximation is constructed as follows. The lagrangian is made only from 
the light fields, which are the six leptons, the five quarks, the photon and the 
gluon. At leading order in the inverse heavy scale, for which we take the W 
boson mass Mw, a single effective operator is added, giving the lagrangian (in 
the so-called charge conserving form of the Fermi operator) 
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C cS = CU") + £qed(«°, ml m° q , 1°, q°, A%) + C QCD (a° s , m% q°, Af) (1) 
+ -y| e°7 Q (l - 7s)/i° x Xyy a (l - 7sK, 

where the superscript denotes bare quantities. The theory is finite after mass 
and coupling (a and a s ) renormalization to all orders in a and a s , and leading 
order in the Fermi constant Gf ~ l/^u^ which is why this parameter is not 
renormalized. 



The matching procedure in the present case consists in requiring that the 
amputated renormalized Green functions 1 of the effective theory be equal 
to the amputated renormalized Green functions of the Standard Model up to 
terms of order C(l/M^) and given order in a and a s 

G S M = G cfi + 0(l/M^), (2) 



which makes the muon decay amplitude the same in both models up to the 
specified order. The Fermi constant is then given as an expansion in a and ct s 

i=0 V £o\v m W 



with Gp = na / (y/2s^ v M^ v ) being the Born level prediction. The quantity 
Ar is customarily used to parametrize the higher order contributions. At the 
one loop level, the matching equation is schematically depicted in Fig. 1. We 
introduced there the decoupling coefficients [13], Z* which are different from 
one in the MS scheme for example, but can be neglected in the on-shell scheme. 
The renormalization constant of the Fermi operator Zq f , although trivial (i.e. 
equal to one), has been included for generality. 

The matching equation, Eq. 2, can be solved in different ways. The apparently 
simplest is to put all light masses and external momenta to zero, and renor- 
malize the wave functions in the on-shell scheme. The right hand side in Fig. 1 
will then consist of only one term, proportional to Gp\ if we use dimensional 
regularization, whereas the left will only have vacuum diagrams with heavy 
masses. Obviously, this situation will persist to all orders. The price to pay 
for this simplicity is the problem of infrared divergent box diagrams, where a 
product of gamma matrices occurs which does not have the form of the Fermi 
operator. In [2,3,4], this product has been defined through Fierz symmetry 
with respect to the last line in the string, which has been implemented, for 

1 This choice is somewhat arbitrary, since one might just as well use full Green 
functions, or Green functions which are one particle irreducible with respect to the 
light fields. 
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Fig. 1. Matching equation at the one loop level. The wavy lines on the left hand 
side represent the three gauge bosons, 7, W and Z, whereas on the right, only the 
photon. The rest of the notation is explained in the text. 




Fig. 2. A three loop diagram, which Fig. 3. A triangular femion loop requir- 
cannot be defined by Fierz symmetry ing special treatment of the 75 matrix, 
with respect to the last line. 

practical reasons, by means of a suitable projection operator. Although suf- 
ficient at the two loop level, this symmetry will not suffice at the three loop 
level, see for example Fig. 2, where due to crossings, there is no last line in this 
sense. It is also not trivial that this procedure is correct even at the two loop 
level. For many topologies, e.g. those that contain a self energy insertion on the 
gauge boson line, one can convince oneself that this is indeed the case, others 
like the nonplanar double boxes in the purely bosonic contributions are not 
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that easy. A highly nontrivial test of the calculation would consist in perform- 
ing the matching without generating spurious infrared divergences 2 . The box 
diagrams being ultraviolet finite can then be calculated in four dimensions, 
avoiding completely the problem of ambiguous gamma matrix definitions. 

In this work we performed the matching by keeping a common mass for all 
the light fermions and evaluating the box diagrams in four dimensions. It 
was necessary to calculate both sides of the matching equation and reexpand 
them subsequently in this common mass. The external wave function renor- 
malization constants were not taken in the on-shell scheme, because this would 
introduce the usual on-shell infrared divergence. On the contrary, they were 
evaluated at zero momentum, which, in practice, is equivalent to renormaliza- 
tion in the MS scheme with nonvanishing decoupling coefficients. Moreover, 
it turned out that it is necessary to have a correct W boson wave function 
renormalization constant, since the box diagrams are not gauge invariant by 
themselves, and in the massive case this constant cancels only in combina- 
tion with vertex diagrams. We used a photon mass regulator, but the MS 
renormalization constant would have been just as good. In the end, complete 
agreement was found with the calculation performed with massless fermions 
and with the projector conserving Fierz symmetry with respect to the last line 
from [2,3,4]. 



3 Results 

A detailed presentation of the methods used to evaluate the bosonic contri- 
butions to Ar( a ) can be found in [4]. The fermionic contributions introduce 
two additional problems. First, some of the two loop vertex diagrams contain 
closed triangular fermion loops as shown in Fig 3. The 75 matrix that occurs 
in the trace has to be correctly defined. We chose the naive dimensional reg- 
ularization scheme [16], with an anticommuting 75 and the four dimensional 
value of the trace of four gamma matrices and 75 

T^tVtVts) = 4ie a ^ 5 . (4) 

This choice is justified by the fact that the nonvanishing contribution of the 
purely four dimensional e tensors is finite. Moreover, it has been checked in [1] 
that the use of the consistent definition of 't Hooft and Veltman [17] gives the 
same result after correction of the Green functions by suitable finite countert- 
erms restoring the Slavnov-Taylor identities. Second, the inclusion of fermions 

2 Note that one could also introduce evanescent operators [14] and perform the 
calculation with vanishing fermion masses and without projection, as for example 
in [15]. This would be a second independent test. 
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results in unstable gauge bosons, which makes a proper definition of their 
masses necessary if gauge invariance of Gp is to be maintained [18]. We use 
the pole mass scheme, where the inverse propagator matrix 

(s-M 2 )5 tJ -U TlJ (s), i,j = W, 1 ,Z, (5) 



is singular in the complex s plane at points which can be parametrized as 
s P = Mp — iMpTp, (6) 



where M P is the mass and T P is the width of the boson. This generates a fixed 
width Breit-Wigner behavior of the total cross section 



as opposed to the running width parametrization actually used by the ex- 
perimental collaborations for the masses and widths of the W and Z bosons 
[19] 

a{s) ~ (s-Mi P y + s ^ cxp /Mi P (8) 



We translate back and forth between the two definition with the help of the 
following relations 

/ F 2 \-l/2 / F 2 \-l/2 



As in [1], we take Y z as experimentally measured, whereas we assume IV to 
be given by the one loop QCD corrected value 

2y/2n I 3tt / 



The complete result for Ar at order a 2 and the partial contributions are 
given in Fig. 4. The top quark mass and the running of the fine structure 
constant are taken from Tab. 1, whereas the masses of the gauge bosons are 
translated from the experimental values given there to the pole mass scheme 
values with the help of Eq. 9, which in this case gives M w = 80.424 GeV and 
M z = 91.1535 GeV. 
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Fig. 4. Complete two loop electroweak contributions to Ar (solid line) together with 
partial corrections: bosonic (dotted line), fermionic (dashed line), light fermionic 
without b quark, but with running of the fine sturcture constant (dash-dotted line) 
and top-bottom (long dashed line). 



input parameter 


value 


source 


M w 


80.451(33) GeV 


[5] 


M z 


91.1876 GeV 


[5] 


m t 


174.3(51) GeV 


[5] 


m b 


4.7 GeV 


[1] 


G, 


1.16637 x 10" 5 GeV" 2 


[20] 


a- 1 


137.03599976 


[5] 


Aa 


0.059228(209) 


[21] 


a s {M z ) 


0.119 


[5] 


r z 


2.4952 GeV 


[5] 



Table 1 

Input parameters with experimental errors, where necessary for the present work. 
The value of m^, is the same as in [1] for comparison purposes. 

In order to compare our result for the fermionic contributions with [1], we 
evaluate the W boson mass from the formula 



Mi 



w 



. 1 1 7T« 
My - + 

'\ 2 \ 4 V2G F M 2 Z 



(1 + Ar), 
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A r ( Q ) + Ar!° 2) 

T ^1 f crm 


_|_ /\ r (a«,) _|_ / \ r («Os) 








i\/T n i 

My/ |1J 


My/ 


A A/f 
ZAM14/ 


M w 


AM W 


[GeVJ 


[GeVj 




[ivie v j 


[GeV] 


[MeV] 


100 


80.3771 


80 3758 


-1.3 


80.3747 


-2.4 


200 


80.3338 


80.3326 


-1.2 


80.3321 


-1.7 


600 


80.2521 


80.2509 


-1.2 


80.2508 


-1.3 


1000 


80.2135 


80.2122 


-1.3 


80.2129 


-0.6 



Table 2 



Comparison of the My/ prediction (third and fifth column) against [1] (second 
column). AMy/ is, in both cases, the shift with respect to the fitting formula. 

with 

Ar = Ar {a) + Ar (aQs) + Ar {aa ^ + a££- ( 12 ) 

We keep a finite b quark mass in Ar^ and Ar' ms ' and take the result for 
Ar' aas ' from [22]. Note also that we do not resum the running of the fine 

structure constant, i.e. Ar^ contains the term +Ao and Arf^ includes 
+Aa 2 . The result is summarized in Tab. 2 for different Higgs boson masses 
from the range from 100 GeV to 1 TeV. We observe a discrepancy of around 
— 1.3 MeV with respect to [1], which comes solely from the differing fermionic 
contributions 3 . 

Inclusion of the bosonic corrections generates an additional variable shift al- 
ready given in [2,3]. As a result, our complete contributions induce a change 
of the My/ prediction by —2.4 MeV for a Higgs boson mass as low as 100 GeV 
(see Tab. 2). Since the bosonic part becomes negative for a heavier Higgs 
boson, this shift reaches -0.6 MeV for M H = 1 TeV. 

In Tab. 3, we include also the recent partial results at three loop order [11], 
i.e. we use 

Ar = Ar(°» + Ar^ + Ar (aa ^ + A ( " 2) - % (^pf ] + ^pf^) ■ (13) 

Together with errors coming from the top quark mass and the running of 
the fine structure constant but without a theoretical error estimate, the M w 

3 The authors of [1] traced a problem in their calculation and after corrections agree 
with our results both for the fermionic and for the Higgs boson mass dependence of 
the bosonic two loop contributions. We checked that all of the remaining corrections 
are the same to required numerical accuracy. 
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C w/ S w n -Pt 




Mw 


M w 


AM W 


M w 


AM W 


[GeVJ 


[GeV] 


[GeV] 


[MeV] 


[GeV] 


[MeV] 


100 


80.3747 


80.375 


0.3 


80.3771 


2.4 


200 


80.3321 


80.3322 


0.1 


80.3358 


3.7 


600 


80.2508 


80.2510 


0.2 


80.2579 


7.1 


1000 


80.2129 


80.2146 


1.7 


80.2231 


10.2 



Table 3 

Additional shift of Mw with respect to the complete prediction from Tab. 2 due 
to inclusion of partial results at order a 3 and a 2 a s from [11]. 
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Fig. 5. The theoretical prediction for the W boson mass, with error from the 
uncertainty of the top quark mass and the running of the fine structure constant, 
against the current experimental value, M^ p . 



prediction is shown against the current experimental result in Fig. 5. 



4 Conclusions 



We have presented a new result for the complete electroweak contributions to 
the lifetime of the muon, which induces a shift in the W boson mass prediction 
as large as —2.4 MeV for a light Higgs boson, of which —1.3 MeV come from a 
discrepancy with the previous calculation of the fermionic contributions [1] and 
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the rest from the bosonic part. The authors of [1] corrected their evaluation 4 
and are now in full agreement with this work. Together with recent results 
at the three loop level [11], this calls for an updated fitting formula. Such a 
formula will be given in a subsequent publication [24]. 
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